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Abstract 

I evaluated the use of global remote sensing techniques for estimating plant leaf 
chlorophyll a + b (C a t,; /(g cm -2 ) and water (C w ; mg cm -2 ) concentrations as 
well as the ratio of C w /C a b with the PROSAIL model under possible distribu- 
tions for leaf and soil spectra, leaf area index (LAI), canopy geometric structure, 
and leaf size. First, I estimated LAI from the normalized difference vegetation 
index. I found that, at LAI values <2, C a b> C w , and C w /C a b could not be reliably 
estimated. At LAI values >2, C a b and C w could be estimated for only restricted 
ranges of the canopy structure; however, the ratio of C w /C ab could be reliably 
estimated for a variety of possible canopy structures with coefficients of deter- 
mination (R 2 ) ranging from 0.56 to 0.90. The remote estimation of the C w /C a b 
ratio from satellites offers information on plant condition at a global scale. 
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Introduction 

Leaf chlorophyll a + b (C a b; fig cm -2 ), dry matter (C m ; 
mg cm -2 ), and water (C w ; mg cm -2 ) concentrations 
provide information on plant physiological status and 
ecosystem functioning (e.g., terrestrial heat, water, and 
C0 2 balances) (Zarco-Tejada et al. 2003). Nevertheless, 
global remote estimates of vegetation status mainly focus 
on leaf area index (LAI; Garrigues et al. 2008), fraction of 
absorbed photosynthetically active radiation (FAPAR; 
Running et al. 2004), phenology (Zhang et al. 2006), leaf 
clumping (Chen et al. 2005), and vegetation height (Sim- 
ard et al. 2011). 



Site-specific studies have clarified the relationships 
between plant canopy spectral reflectance and C a b (Zarco- 
Tejada et al. 2004; Gitelson et al. 2005; Darvishzadeh et al. 
2012; Si et al. 2012), C m , (Fourty and Baret 1997) and C w 
(Bowyer and Danson 2004; De Santis et al. 2006; Zarco-Tej- 
ada et al. 2003). However, global or regional relationships 
between these factors are not clearly understood because 
other characteristics of the plant canopy impact those rela- 
tionships. In this note, I evaluated the use of global remote 
sensing techniques for estimating C ab , C m , C w , and the CJ 
C a b ratio using the PROSAIL model (Jacquemoud et al. 
2009) at varying possible distributions of leaf and soil 
spectra, LAI, canopy geometric structure, and leaf size. 



2544 



© 2012 The Author. Published by Blackwell Publishing Ltd. This is an open access article under the terms of the Creative 
Commons Attribution Non-Commercial License, which permits use, distribution and reproduction in any medium, provided 
the original work is properly cited and is not used for commercial purposes. 



K. Kushida 



Remote Sensing Leaf Water-Chlorophyll Ratio 



Table 1. Input parameters for the PROSAIL model. 



Unit 



Value/function 



ln(C ab ) 
ln(CJ 
ln(C w ) 

The coefficient of correlation 

between ln(C a b) and ln(C m ) 
The coefficient of correlation 

between ln(C w ) and ln(C ab ) 

or ln(CJ 
The parameter characterizing 

the leaf mesophyll 

structure, W 
Soil reflectance 

Clumping index, fl 

Leaf angle distribution, LAD 



Canopy hotspot 

parameter, 5\ 
Leaf area index, 

LAI = LAI r • n 
Solar incident zenith angle, 0 S 
View zenith angle, 9 V 
Solar illumination specular 

ratio, r sd 



Error function of the 
atmospheric 
correction 



ng cm- 2 W(3.79, 0.35 2 )* 
mg cm" 2 N(1 .57, 0.42 2 )* 
mg cm" 2 W(2.32, 0.47 2 )* 
0.56 



N(1-7, 0.2 2 )* 



JHU-SL based 

statistic model 
Any (fixed at 0.7 for 

LAI calculation) 
Erectophile, spherical, 

plagiophile, uniform, 

extremophile, and 

planophile 
0.0001, 0.001, 0.01, 

and 0.1 
0-10 (10~ 6 interval) 

25 and 50 
0 

Constant (0.81, 0.91, 
0.95, 0.98, and 1.0 
in MODIS bands 
3, 4, 1 , 2, and 5-7) 

W(0, (0.005 + 0.05p) 2 )* 



*N(ji,a ) denotes the normal distribution with mean \i and variance a . 

Materials and Methods 

I used the PROSAIL model, which is a combination of the 
SAIL (Verhoef 1984) and PROSPECT (Jacquemoud and Ba- 
ret 1990) models, for calculating the relationships between 
the top-of-the-atmosphere (TOA) canopy spectral reflec- 
tance and C a b, C m , C w , and C w /C a b- The input parameters of 
the PROSAIL model are listed in Table 1, and I noted TOA 
canopy spectral reflectance from the model outputs. 

I calculated the relationship between the normalized dif- 
ference vegetation index (NDVI) and LAI (m 2 m~ 2 ) follow- 
ing Kushida and Yoshino (2010). I used the TOA canopy 
spectral reflectance values at the Moderate Resolution Imag- 
ing Spectroradiometer (MODIS) red (620-670 nm; R R (%)) 
and near infrared (841-876 nm; ^ N i R (%)) bands to calculate 
NDVI as follows: 



NDVI 



Rnir + Rr 



(1) 



For the estimated LAI value ranges of 0.95-1.05, 1.95- 
2.05, 2.95-3.05, 3.95-4.05, 4.95-5.05, 5.95-6.05, and 6.95 
-7.05, I calculated the relationships between canopy spec- 
tral reflectance and C a b, C m , C w , and C w /C a b- 

I assumed that C a b, C m , and C w had lognormal distri- 
butions based on the leaf optical properties experiment 
(LOPEX93), which provides values for leaf pigment 
and water content of 70 leaf samples that represent 
approximately 50 species of woody and herbaceous plants 
(Hosgood et al. 1995). The means and standard 
deviations (SD) of the lognormal distributions and the 
correlations and ranges of the variables were also deter- 
mined using the data values in LOPEX93. All calculations 
were carried out under variable soil spectrum conditions. 
The soil reflectance distribution for each of the bands was 
assumed to be lognormal or normal, based on the Johns 
Hopkins University Spectral Library (JHU-SL; Baldridge 
et al. 2009). The means and SD of the distributions and 
the correlations and ranges of the variables were deter- 
mined using the data values in JHU-SL. 

The clumping index (Q; Chen et al. 2005) was incor- 
porated in the model by setting the parameter at 0.7 to 
express the leaf clumping effect. Using this parameter, I 
adjusted the LAI value such that, in a plant canopy with 
an initial LAI value of LAI T and an f2 value of fl, the 
adjusted LAI value became LAI T ■ fllQ.7. The leaf angle 
distribution (LAD) was fixed as erectophile, spherical, 
plagiophile, uniform, extremophile, and planophile. The 
canopy hotspot parameter (Si), which is equal to the ratio 
of the correlation length of leaf projections in the 
horizontal plane and the canopy height, was fixed at 
0.0001, 0.001, 0.01, and 0.1. I used these ranges of LAD 
and Si values to represent global distributions in these 
parameters. The solar incident zenith angle (9 S ) was 
fixed at 25° for all model evaluations except when model 
sensitivity to this parameter was evaluated, and, in this 
case, 9 S was increased to 50°. The specular ratio to the 
total solar illumination (r sd ) was set at constant values 
representing typical atmospheric conditions on a clear 
day for each of the bands. I assumed that the error 
function of the atmospheric correction had an indepen- 
dent normal distribution with a mean of 0 and SD of 
0.005 + 0.05p (no unit reflectance), where p is the 
reflectance value at a given spectral band (Vermote and 
Kotchenova 2008). 

In the calculation of each of combinations of the LAD 
and Si types, LAI values were provided from 0 to 10 at 
10~ 6 intervals. I used pseudo-random numbers to express 
the lognormal and normal distributions of the leaf spec- 
tral parameters and soil reflectance. I calculated 10 7 cases 
to obtain the relationships between C a b, C m , and C w and 
their associated spectral reflectances. 
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I used the TOA canopy spectral reflectance values at the 
MODIS green (545-565 nm; R G (%)), near infrared (841- 
876 nm; -Rnir(°/°))> anQ shortwave infrared (1628- 
1652 nm; -Ri 6 4o(%)) bands to estimate C a b, C m , and C w , 
respectively. The R G , -R N i R , and -Ri 6 4o bands correspond to 
absorption bands of chlorophyll a + b, dry matter, and 
water and dry matter combined, respectively. Absorptions 
of carotenoids and anthocyanins also concern with _R G ; 
however, in general, most of the leaf absorption at the 
green band corresponds to chlorophyll a + b. That was 
because the carotenoids concentration and C a b have a high 
positive correlation for a variety of plant species and the 
anthocyanins concentration appears when the leaves of 
deciduous trees turned red in autumn. For the estimation 
of C w , the estimated C m value from _R N ir was multiplied by 
the ratio of the specific absorption coefficient of water 
(cm 2 mg) to that of dry matter (cm 2 mg), which is 0.787, 
and then removed this value (/i(-Rnir)) from the estimated 
C w because the PROSPECT model shows that both C w and 
C m contribute to i?i 6 4o- For the estimation of C w /C a b, the 
estimated C w value was divided by the estimated C a b value. 

I divided the ranges of the values of the TOA spec- 
tral reflectance or the abovementioned spectral indices 
into 8-12, and then calculated the average and SD of C a b, 
C m , C w , and C w /C a b for each of the divided ranges. To 
express the relationships in mathematical formulas for 
each of the estimated LAI value ranges, I used regression 
equations in the form: 

y = ax~\ (2) 

where a and b are constant values, x is the spectral reflec- 
tance (_R G , -Rnir> or R 1640 + /;(^nir))> and y is the leaf 
constituent (C a b, C m , or C w ). For the ratio C w /C a b for 
each of the estimated LAI value ranges, I used the regres- 
sion equation in the form: 

= flli<G KI640 - a 2 Rc «NIR , (3) 

Cab 

where a x , a 2 , b\, b 2 , c u and c 2 were constant values. The 
coefficient of determination (_R 2 ) was defined as: 



where E is the root mean square error (RMSE) from the 
regression and SD is the standard deviation of the sam- 
ples. Standard deviation for C a b, C m , C w , and C w /C a b was 
16.4 /.(g cm~ 2 , 2.22 mg cm~ 2 , 5.63 mg cm~ 2 , and 0.174, 
respectively. 

Results 

The LAI was estimated from the NDVI with R 2 of 0.39 for 
all combinations of the LAD and Si types together when 




NDVI 

Figure 1. Relationship between normalized difference vegetation 
index and estimated leaf area index for all combinations of leaf angle 
distribution and S| types (8 S = 25°). The solid and dotted curves 
denote the average and the average ±SD, respectively, for each 
estimate. 

6 S = 25° (Fig. 1). The SD of the estimated LAI values over 
the intervals 0.95-1.05, 1.95-2.05, 2.95-3.05, 3.95-4.05, 
4.95-5.05, 5.95-6.05, and 6.95-7.05 were 0.4, 0.9, 1.6, 2.1, 
2.2, 2.1, and 1.8, respectively. Similarly, the LAI was esti- 
mated form the NDVI with R 2 of 0.30 when 8 S = 50°. 

I found a weak relationship between C m and -R N ir 
and a strong relationship between C a b, C w , and C w /C a b 
and their associated canopy spectral reflectances at the 
LAI values >2 when all the LAD and Si types were 
equally probable in one pixel of a remotely sensed image 
(Fig. 2). At the LAI values of 1, the R 2 of C ab , C w , and 
C w /C a b estimations were <0.22. The LAI estimates of 2, 
4, and 6 in Figure 2 correspond to the estimated LAI 
ranges of 1.95-2.05, 3.95-4.05, and 5.95-6.05, respec- 
tively. 

However, for each of the estimated LAI ranges, as the 
LAD became more vertical or as Si decreased, C a b, C m , 
and C w decreased for the same canopy spectral reflectance 
values (Fig. 3). As the relationships between leaf constitu- 
ents and associated canopy reflectances were dependent 
on the LAD and Si types, it was difficult to estimate C a b 
and C w from the canopy spectral reflectances when 
different LAD and Si types existed in the focal region of 
analysis. R p 2 and R 2 in Figure 3 were the coefficients of 
determination under erectophile and Si = 0.0001 parame- 
terization and under planophile and Si = 0.1 parameteri- 
zation, respectively. In contrast, for the estimation of the 
C w /C a b ratio, estimation equations were independent of 
the LAD and Si types (Fig. 3). The relationships between 
the estimated and observed C w /C a b ratios for all combina- 
tions of the LAD and Si types (Table 1) fell within the 
two regression curves for erectophile and Si = 0.0001 
parameterization and for planophile and Si = 0.1 parame- 
terization (Fig. 3). 
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1 

25 

R 1640 +f i ( R NIr)( 0/, °) 



50 0 



— 

■a 0.5- 




R 2 = 0.58 



— V 

0.5 




R 2 = 0.72 
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R 2 = 0.17 
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Figure 2. fi G versus C a b, /?nir versus C m , /?i64o + Wnir) versus C w , and estimated C„/C a b versus observed Cv/Cab for all combinations of leaf angle 
distribution and S\ types and for estimated leaf area indices of 2, 4, and 6 (8 S = 25°). The solid and dotted curves denote the average and the 
average ±SD, respectively, for each estimate. 



I found that the models were not sensitive to the 
value used for 6 S . For an LAI value of 1, the R 2 of C a b, 
C w , and C w /C a b estimates were <0.30, and the relation- 
ship between C m and the canopy spectral reflectance was 
weak. I could not estimate C a b and C w when different 
LAD and Si types existed in the focal region of analysis 
because the equations used to estimate C a b and C w were 
dependent on the LAD and Si types. However, the esti- 
mation equations for C w /C a b were independent of the 



LAD and Si types. Coefficients of the regression equa- 
tions of C w /Qb for all combinations of LAD and Sj 
types together in the form of equation (3) and the R 2 
are shown in Table 2. 

Discussions and Conclusion 

I evaluated global estimates of C a b, C m , C w , and C w /C a b 
from TOA broadband spectral reflectance using the PRO- 
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Figure 3. R G versus C a b, Rnir versus C m , /?i64o + Wmr) versus C w , and estimated C„/C a b versus observed C„/C a t, for estimated leaf area indices of 
2, 4, and 6 (6 S = 25°). Lines denote the averages of all combinations of leaf angle distribution (LAD) and S\ types (solid lines), erectophile and 
S| = 0.0001 (dotted lines), and planophile and S, = 0.1 (dashed lines). R 2 e represents R 2 when the LAD is erectophile and 5, = 0.0001, and R 2 P 
represents R 2 when LAD is planophile and Si = 0.1. 



SAIL model for possible distributions of leaf and soil 
spectra, LAI, canopy geometric structure, and leaf size. 
For LAI values <2, C a b, C w , and C w /C a b had weak rela- 
tionships with their associated spectral reflectances. For 



LAI values greater than 2, C a b and C w could be reliably 
estimated only for certain ranges of LAD and Si types 
while the ratio of C w /C a b could be reliably estimated for 
all possible canopy structures with an R 2 ranging from 
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Table 2. Coefficients of the equations for estimating CJC^. 



8 S 


LAI 


91 




b 2 


^2 


Ci 


c 2 


R 2 


25 


2 


12.231 


1.003 


1.659 


0.177 


0.859 


0.691 


0.58 




3 


8.660 


0.949 


1.541 


0.349 


0.858 


0.847 


0.67 




4 


6.632 


0.903 


1.446 


0.443 


0.781 


0.845 


0.72 




5 


5.392 


0.866 


1.373 


0.665 


0.682 


0.876 


0.77 






4 604 


0 841 


1 326 


1 300 


0 587 


0 981 


0 83 




7 


4.179 


0.784 


1.286 


3.604 


0.529 


1.208 


0.89 


50 


2 


12.120 


0.978 


1.645 


0.239 


0.898 


0.785 


0.61 




3 


8.776 


0.914 


1.525 


0.400 


0.881 


0.888 


0.67 




4 


6.745 


0.870 


1.430 


0.532 


0.807 


0.900 


0.73 




5 


5.727 


0.831 


1.369 


0.937 


0.706 


0.971 


0.78 




6 


4.932 


0.810 


1.327 


2.136 


0.614 


1.117 


0.84 




7 


4.549 


0.726 


1.281 


7.318 


0.535 


1.391 


0.90 



0.56 to 0.90. The estimation equation of C w /C a b from the 
spectral reflectance and the R 2 value were dependent on 
the LAI values, but independent of the LAD and S[ types. 

Levels of C a b, C m , and C w can be indicative of leaf 
physiology and plant condition, and attempts have been 
made to estimate these values with remote sensing 
applications (Ustin et al. 2009). In previous site-specific 
studies, C a b, C m , and C w were successfully estimated 
using this methodology (Fourty and Baret 1997; Bowyer 
and Danson 2004; Zarco-Tejada et al. 2004, 2003; Gitel- 
son et al. 2005; De Santis et al. 2006; Darvishzadeh et al. 
2012; Si et al. 2012); however, the use of remote sensing 
to estimate these values at regional and global scales has 
not been reported. My research suggests that estimates 
of C a b, C m , and C w are dependent on LAD and Si types, 
and these types often vary across regional to global 
scales. As techniques for estimating LAD and Si types 
through remote sensing have not been established, a 
generalized estimation of C a b, C m , and C w across broad 
spatial scales is difficult except in cases where specific 
LAD and Si types can be inferred. In contrast, this study 
shows that the estimation of the C w /C a b ratio through 
remote sensing techniques is generally possible across 
regional and global scales. The estimation of this ratio 
through remote sensing has not been previously consid- 
ered as an indicator of plant canopy condition because 
the physiological meaning of the C w /C a b ratio has been 
less well studied than those of C a b, C m , and C w . However, 
I found that the C w /C a b ratio had a stronger relationship 
with its associated spectral band than C a b, C m , and C w 
had to their associated bands. 

Previous studies described the meaning behind and 
variations in C a b, C w , and C w /C a b- The value of C w /C a b is 
specific to plant species, although, in general, the ratio 
slightly decreases with spring sprout and increases with 
autumn defoliation (Gond et al. 1999; Ceccato et al. 
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2001). For example, Scots pine, lodgepole pine, sun flow- 
ers, and sugar beets have high C w /C ab values (1.0-1.8); 
poplars, oaks, and rhododendrons have moderate ratio 
values (0.5-0.7); and maize and rice have low values (0.1- 
0.2; Ceccato et al. 2001; Gond et al. 1999; Hosgood et al. 
1995). The value of C w /C a b also generally increases when 
a plant responds to stressors related to water deprivation 
(Zhang and Kirkham 1996; Guerfel et al. 2009), heat 
(Jeon et al. 2006), chilling (Bacci et al. 1996; Jeon et al. 
2006; Korkmaz et al. 2010), high light conditions (Jagtap 
et al. 1998), ultraviolet rays (Alexieva et al. 2001), high 
salinity (Jaleel et al. 2008; Dogan 2011), and heavy-metal 
contaminants (Anuradha and Rao 2009). The increase in 
the ratio as a result of plant stress is caused by a greater 
decrease in leaf chlorophyll compared with leaf water. 

Therefore, a change in the C w /C a b ratio through time is 
the result of either changes in species composition or 
changes in the response of plants to stress. The former gen- 
erally occurs at a yearly to decadal scale, whereas the latter 
occurs at a daily to monthly scale. Although disturbances 
such as wildfires, insect attacks, and deforestation can 
cause an immediate change in species composition, an 
analysis of changes in LAI and NDVI can help to distin- 
guish between the causes of changes in C w /C a b- Thus, the 
remote estimation of the C w /C a b ratio from satellites offers 
information on plant status at a global perspective. 

Other than leaf chlorophyll a + b, leaf anthocyanins 
absorb the green light and generate the errors in the esti- 
mation of the C w /C a b ratio. Leaf anthocyanins appear 
when the leaves of deciduous trees turned red in autumn 
or under chilling stress (Bacci et al. 1996). When a 
change in the estimated C w /C a b ratio possibly appeared in 
such cases, an analysis of changes in the estimated LAI 
and NDVI can help to distinguish between the causes of 
changes in the estimated C w /C ab , as the autumn colora- 
tion reduces NDVI (Zhang et al. 2012). 

This study modeled the global vegetation as to obey 
PROSAIL model with parameters shown in Table 1. The 
R 2 of the LAI estimation from NDVI with globally fixed 
equations were 0.30-0.39, whereas the .R 2 of the C w /C a b 
ratio estimation in the way I presented in this note was 
0.56-0.90, for the estimated LAI values >2. This indicates 
that the global remote estimation of the C w /C a b ratio is 
more reliable than that of LAI. 
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